Astroparticle Physics 45 (2013) 56-70 



ELSEVIER 


Contents lists available at SciVerse ScienceDirect 

Astroparticle Physics 

journal homepage: www.elsevier.com/locate/astropart 



The Astrophysical Multimessenger Observatory Network (AMON) 



CrossMark 


M.W.E. Smith a,b ’*, D.B. Fox c,b , D.F. Cowen a,c,b , P. Meszaros a,c,b , G. Tesic a,b , J. Fixelle c , I. Bartos d , 

P. Sommers a,c,b , Abhay Ashtekar a b , G. Jogesh Babu e,c , S.D. Barthelmy f , S. Coutu a,b , T. DeYoung a,b , 
A.D. Falcone c,b , Shan Gao a,b , B. Hashemi 3 , A. Homeier g , S. Marka d , B.J. Owen a,b , I. Taboada 11 


a Department of Physics, Pennsylvania State University, University Park, PA 1 6 802, USA 
b Institute for Gravitation and the Cosmos, Pennsylvania State University, University Park, PA 1 6802, USA 
c Department of Astronomy & Astrophysics, Pennsylvania State University, University Park, PA 16802, USA 
6 Department of Physics & Columbia Astrophysics Laboratory, Columbia University, New York, NY 10027, USA 
e Department of Statistics, Pennsylvania State University, University Park, PA 1 6802, USA 
'NASA Goddard Space Flight Center, Code 661, Greenbelt, MD 20771, USA 
g Physikalisches Institut, Universitat Bonn, Nussallee 12, D-53115 Bonn, Germany 

h Center for Relativistic Astrophysics, School of Physics, Georgia Institute of Technology, Atlanta, GA 30332, USA 


ARTICLE INFO 


ABSTRACT 


Article history: 

Received 3 December 2012 

Received in revised form 1 February 2013 

Accepted 15 March 2013 

Available online 28 March 2013 


Keywords: 

High-energy astrophysics 
Gravitational radiation 
Neutrinos 
Cosmic rays 
Gamma-ray bursts 
Supernovae 


We summarize the science opportunity, design elements, current and projected partner observatories, 
and anticipated science returns of the Astrophysical Multimessenger Observatory Network (AMON). 
AMON will link multiple current and future high-energy, multimessenger, and follow-up observatories 
together into a single network, enabling near real-time coincidence searches for multimessenger astro- 
physical transients and their electromagnetic counterparts. Candidate and high-confidence multimessen- 
ger transient events will be identified, characterized, and distributed as AMON alerts within the network 
and to interested external observers, leading to follow-up observations across the electromagnetic spec- 
trum. In this way, AMON aims to evoke the discovery of multimessenger transients from within obser- 
vatory subthreshold data streams and facilitate the exploitation of these transients for purposes of 
astronomy and fundamental physics. As a central hub of global multimessenger science, AMON will also 
enable cross-collaboration analyses of archival datasets in search of rare or exotic astrophysical 
phenomena. 

© 2013 Elsevier B.V. All rights reserved. 


1. Introduction 

We stand at the dawn of multimessenger astrophysics - a quest 
to use the messenger particles of all four of nature’s fundamental 
forces to explore the most violent phenomena in the universe. 
Observatories first imagined a generation ago are finally being 
realized, including the Advanced L1G0 [1] and Virgo [2] gravita- 
tional-wave detectors, the ANTARES [3] and IceCube [4] high- 
energy neutrino observatories, and the Pierre Auger Cosmic Ray 
Observatory [5]. On the ground and in space, they are comple- 
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mented by high-energy observatories including the Swift [6] and 
Fermi [7] satellites, the HESS [8], VERITAS [9], and MAGIC [10] 
TeV gamma-ray telescopes, and the HAWC [11] TeV gamma-ray 
observatory. 

Collectively, these facilities promise the first detections of grav- 
itational waves and high-energy cosmic neutrinos, the resolution 
of the mystery surrounding the origins of ultrahigh-energy cosmic 
rays, and a new window into the formation and evolution of black 
holes. Given their nature as first-generation facilities, however, the 
sensitivities of the non-electromagnetic observatories are naturally 
limited, with rates of detection for transient events of publishable 
significance known or expected to be low, perhaps a handful per 
year (or, in the case of Advanced LIGO and Virgo at design sensitiv- 
ity, a few dozen [1]). 

During the intervals prior to and between detection of these 
rare high-significance events, the multimessenger facilities will 
be buffeted by signals from a far greater number of lower-signifi- 
cance events that will be statistically indistinguishable from back- 
ground or noise processes. Such “subthreshold” events are, by 
definition, unrecoverable as astrophysical signals in the data 
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stream of any individual facility. However, if they are accompanied 
by a subthreshold signal in another multimessenger channel they 
can be identified, and potentially achieve high significance, via 
careful coincidence analysis of the data streams from multiple 
facilities. 

In this paper we present the scientific case for the Astrophysical 
Multimessenger Observatory Network and describe its important 
elements. AMON will weave together existing and forthcoming 
high-energy astrophysical observatories into a single virtual sys- 
tem, capable of sifting through the various data streams in near 
real-time, identifying candidate and high-significance multimes- 
senger transient events, and providing alerts to interested 
observers. 

As we show, AMON will enable a significant enhancement in 
the effective aggregate sensitivity of the world’s leading multimes- 
senger facilities for a small fraction of the facilities’ total cost, pro- 
vide the first near real-time alerts for multimessenger transient 
sources, and simplify the mechanics and politics of cross-collabo- 
ration analyses for all partners. As such, we believe AMON repre- 
sents a natural next step in the extension of the global 
astronomical community’s vision beyond the electromagnetic 
(EM) spectrum. 

The development of AMON is currently underway. Signatories 
to the AMON Memorandum of Understanding 1 (MOU) include the 
IceCube and ANTARES neutrino observatories, the VERITAS and 
HAWC TeV y-ray observatories, and the Swift orbital telescope. 
Exploratory discussions with the LIGO (including GEO-600), and Fer- 
mi scientific collaborations have led to signed letters of commit- 
ment, with negotiations toward MOU signatures from all parties 
ongoing. Discussions have also been initiated with candidate fol- 
low-up facilities including ROTSE-III [12] and the Palomar Transient 
Factory [13], with the goal of bringing these observatories into the 
collaboration prior to or shortly after the commencement of real- 
time AMON alert operations. 

AMON is structured as an open and extensible network, with an 
MOU that allows straightforward incorporation of new triggering 
and follow-up facilities. Collaborations interested in the scientific 
goals of AMON, and with useful triggering or follow-up capabilities 
to contribute, are encouraged to contact the authors for informa- 
tion about joining AMON. First versions of the necessaiy support- 
ing hardware and software infrastructure for AMON are being 
installed at Penn State, and initial analyses on archival and simu- 
lated real-time data streams will get under way shortly thereafter, 
as a means of preparing to bring the first set of triggering facilities 
on-line within the next year. 

The paper is structured as follows: In Section 2 we provide an 
overview of the scientific opportunity for AMON which motivates 
our efforts. In Section 3 we discuss the elements of AMON, includ- 
ing the technical and operational protocols that we propose to 
adopt in linking the partner facilities, and the algorithms that we 
will use to identify coincident signals. Section 4 presents detailed 
simulations of multimessenger transient sources, and additional 
theoretical case studies, which illustrate the gains that stand to 
be realized by AMON. Section 5 presents our summary and 
conclusions. 


2. The AMON science opportunity 

AMON is intended to contribute in several ways to the first dec- 
ade of multimessenger astronomy. A common focus of these ap- 
proaches is on multimessenger transient events that are 
observed as coincident (potentially subthreshold) signals in the 
data of AMON partner facilities corresponding to two or more dis- 
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tinct types of messenger particle. In this section, we briefly review 
the strongest candidates for these multimessenger transient 
sources and the current theoretical expectations for their proper- 
ties, rates, and broader implications for physics and astrophysics. 

2 A. High-luminosity gamma-ray bursts 

As the most violent cataclysms known, and as the sites of the 
most highly-relativistic cosmic outflows, gamma-ray bursts (GRBs) 
have long been considered likely multimessenger transient 
sources. Recent progress in distinguishing the distances and ener- 
getics of GRBs and understanding the likely nature of their progen- 
itors [14] means that we can review the expected multimessenger 
signals of the different varieties of GRB on a case by case basis. 

The traditional or “classical” long-duration, high-luminosity 
gamma-ray bursts (HL-GRBs) are the most common type of GRB 
detected by satellite experiments, being observed as bright sec- 
onds- to minutes-long bursts of y-radiation from high-redshift, 
zi> 1. HL-GRBs are believed to arise when a massive star 
(M > 25 M 0 ) undergoes core collapse to a black hole (BH); confir- 
mation of this “collapsar” model [15] for the HL-GRBs has been 
most dramatically provided by spectroscopic observations of ensu- 
ing type Ibc supernovae (SNe), with other lines of evidence also 
contributing [16]. 

In the collapsar model, formation of a high angular-momentum 
BH and accretion of residual gas through an accretion disk pro- 
duces a relativistic jet. In “successful" GRBs the jet pierces the stel- 
lar envelope, accelerates to high Lorentz factor, and radiates 
gamma-rays for tens of seconds, providing a bright electromag- 
netic trigger for observers within the jet collimation angle. The typ- 
ical HL-GRB jet energy of E ~ 10 51 — 10 52 erg, collimated within an 
angle of ( 0 ) ~ 5°, roughly 1/250 of the sky, yields the observed iso- 
tropic-equivalent energies of E 0 = E(4n/Qj) ~ 10 53 — 10 54 erg [17], 
For a burst duration T~ 10 s, this corresponds to an isotropic- 
equivalent luminosity L y ~ E 0 /T = 10 52 — 10 53 erg s _1 . Given a 
250:1 jet collimation factor, the nearest off-axis bursts are antici- 
pated to be located 250 1/3 ~ 6 times closer than the nearest on- 
axis bursts. Off-axis bursts, defined as bursts whose bright initial 
y-ray emission does not illuminate observers at Earth, may still 
be observed via their less luminous shock-breakout emission, their 
prompt gravitational wave (GW) signal, or their subsequent after- 
glow and/or supernova components. 

In the standard internal-shock model, the gamma-ray emission 
of HL-GRBs is produced by internal shocks at a dissipation radius 
r d ~ r 2 cdt~ (3 x 10 13 cm)^ 5 i5t_2 [18], where )/ 25 = r/300 is 
the bulk Lorentz factor divided by 300 and 5t_ 2 = <5t/l 0 2 s is the 
variability time of the central engine in hundredths of a second. 
These shocks Fermi-accelerate electrons, which produce gamma- 
rays via synchrotron and inverse Compton scattering, and are then 
boosted by the bulk motion of the relativistic outflow. For a discus- 
sion of alternative models see Ref. [19]. 

The same shocks responsible for electron acceleration and y-ray 
emission should also accelerate protons, leading to photo-pro- 
duced pions which in turn yield high-energy cosmic neutrinos 
(CN) and y-rays from charged and neutral pion decays, respectively 
[20,21 ]. Initial assumptions were that the energy in relativistic pro- 
tons would be of the same order as the energy emitted in y-rays, 
E p ~ (1 /f e )Ey, where / e < 1 is the fraction of proton energy given 
to electrons (and observed as GRB photons), so that the optical 
depth T py determines the CN and TeV gamma-ray luminosity, 
Liev ~ ZpyL p . The CN to y-ray flux ratio expected in this model has 
been quantified [22] and used by IceCube [23-25] to set limits 
which are already a factor of five below the naive predictions. 
However, careful consideration of the underlying physics [26-28] 
suggests that IceCube observations will need to continue for sev- 
eral more years before the somewhat reduced CN fluxes of more 
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realistic models are tested. Ultimately, detection of GRB-related CN 
is expected if HL-GRBs contribute a significant fraction of the high- 
est-energy cosmic rays. 

Regardless of the total energy release of the HL-GRBs, which 
may well be an order of magnitude or more above the beaming- 
corrected y- ray energies E < 10 52 erg, no significant GW emission 
will be produced if the core collapse, jet production, and burst pro- 
cesses maintain approximate axisymmetry throughout. However, 
it is possible for a rapidly rotating core and accretion disk to devel- 
op bar and/or fragmentation instabilities which could result in sub- 
stantial GW emission. In the most optimistic case [29] the resulting 
GW signal is periodic and roughly as strong as the signal from neu- 
tron star binary mergers, and hence, visible out to similar distances 
(hundreds of Mpc for the advanced facilities) using ground-based 
GW detectors. 

While some 3D simulations of stellar collapse suggest a frag- 
mentary process, and a correspondingly weaker GW energy release 
Ecw ~ 1CT 7 M p c 2 ~ 10 47 erg [30], these have not included relativis- 
tic effects that are known to contribute to large bar asymmetries 
[31], as manifested in other simulations [32]. When present, these 
instabilities result in GW emission comparable to that found in 
optimistic analytical calculations [29]. 

2.2. Low-luminosity gamma-ray bursts 

Low-luminosity GRBs (LL-GRBs) are underluminous long-dura- 
tion GRBs, having substantially lower isotropic-equivalent ener- 
gies, E 0 ~ 10 49 — 10 50 erg, than the HL-GRBs. Because of their 
lower y-luminosities, they are currently only detected at low red- 
shift, z < 0.5, where they provide the bulk of the observed GRB- 
SN Ibc associations [16]. Estimates of the LL-GRB rate suggest that 
they occur at g 100 times the rate of HL-GRBs while constituting 
<1% the rate of type Ibc SNe overall [33-36]. 

The y-ray emission of LL-GRBs may be due, in some or all cases, 
to a relativistic shock breakout [37,38] and hence may not require 
the Lorentz factor r S 100 jet needed to explain the high-luminos- 
ity, high-variability y-emission of HL-GRBs. Such shock breakout 
emission would also be expected to be uncollimated (uniform over 
the sky) or nearly so. 

Predictions for the CN emissions of LL-GRBs have been explored 
by Murase et al. [39] and Gupta and Zhang [40] for relativistic jet 
models and by Kashiyama et al. [41] for shock breakout models. 
In relativistic jet models, a straightforward approach scales CN 
fluxes from HL-GRB predictions according to their y-ray luminosi- 
ties, and anticipates similar spectra with El<f> v peaking at PeV-EeV 
energies. Given current evidence for less-relativistic jets in LL- 
GRBs, however, these predictions may be optimistic. 

Shock breakout models [41 ], by contrast, predict softer neutrino 
spectra peaking in the TeV-PeV range, with luminosities for refer- 
ence events that would make individual bursts detectable to Ice- 
Cube within a horizon of D < 10 Mpc. The LL-GRB rate within 
this horizon is thought to be small, r~ 0.002 yr _1 , but prospects 
for detection could be improved by stacking analyses using ~doz- 
ens of more distant events identified by their prompt high-energy 
EM signature. 

Similarities in explosion energy, ejecta velocity, and synthe- 
sized nickel mass for the type Ibc supernovae of LL-GRB and HL- 
GRB events suggest that the details of core collapse for the two 
event classes, and their GW emissions, may be similar. However, 
given the critical role that high values of core angular momentum 
are thought to play in powering HL-GRB jets, and the need for 
accretion or disk instabilities to power significant GW emission, 
the relative weakness of LL-GRB jets may suggest correspondingly 
reduced prospects for GW emission. If the GW emissions are com- 
petitive with those of binary neutron star-neutron star (NS-NS) 
mergers in even some cases, this would make the LL-GRBs with 


their relatively isotropic y-ray emissions a highly-promising target 
population for ground-based GW detectors. 

2.3. Short-hard gamma-ray bursts 

Short-hard gamma-ray bursts (SHBs), apart from their shorter 
durations and somewhat harder spectra, are observationally simi- 
lar to HL-GRBs. Their harder spectra may indicate higher bulk Lor- 
entz factors [18], while their reduced durations, t b < 2 s, suggest 
shorter accretion timescales compared to those of the HL-GRB col- 
lapsars. A softer “extended emission” episode lasting for 
30— 100 s after the burst itself is present in about a third of 
SHBs, accounting for 5-20% of the total energy [42]. 

Thanks to afterglow discoveries and host galaxy identifications 
of the Swift era, consensus now holds that SHBs are likely due to 
compact binary (NS-NS or NS-BH) mergers [43,14]. If so, these 
events are associated with progenitor systems that are observed 
within our own Galaxy (as relativistic NS-NS pulsar binaries), have 
merger rates that can be estimated from population synthesis 
modeling (e.g. [44]), and should emit strong GW emission in a 
highly-calculable “chirp” waveform across the frequency range of 
ground-based GW observatories just prior to merger. Indeed, bin- 
ary mergers are the primary extragalactic target population for 
next-generation LIGO and VIRGO, with expected event rates at full 
design sensitivity expected to reach dozens annually [1]. 

Afterglow observations for a handful of Swift-detected SHBs 
show evidence for typical collimation angles (9) > 6°, correspond- 
ing to beaming fractions of 200:1 (similar to that for HL-GRBs) or 
less [45], These beaming corrections adjust the observed isotro- 
pic-equivalent energies of the bursts, E 0 > 10 51 erg to lower in- 
ferred jet energies E ~ 10 48 — 10 5 ° erg. As such the nearest off-axis 
mergers are predicted to be observed at 200 1/3 ~ 6 times closer 
distance than the nearest on-axis events, where they may be de- 
tected as GW inspiral signals without SHB counterparts. 

Because of their very similar y-ray characteristics and inferred 
jet properties, it is assumed that the radiation physics of the HL- 
GRBs and SHBs are similar, apart from possible slight differences 
due to the SHBs’ reduced durations and potentially larger Lorentz 
factors and collimation angles [18]. CN production models may 
thus conservatively be carried forward by assuming similar CN to 
y-ray flux ratios as for the HL-GRBs. 

In terms of GW emission, SHBs in the compact binary merger 
model represent a “dream scenario”, as systems in which the 
GW waveform is known to be strong (detectable to hundreds of 
Mpc by advanced detectors) and calculable, and moreover, has al- 
ready been implemented in data analysis systems via matched-fil- 
ter algorithms that will be run in real-time by LIGO and VIRGO, 
enabling the observatories to distribute GW “inspiral alerts” to 
interested observers. The primary questions for the near future, 
then, are whether the GW inspiral signatures of these mergers will 
be detected, at what rate, and whether or not a coincident y-ray or 
other electromagnetic signature can be identified. Unless our 
understanding is rather radically misplaced, all three questions 
are likely to be resolved once the era of the advanced GW facilities 
is fully under way. 

2.4. Choked jet supernovae 

The choked jet supernova, a theoretical scenario, represents the 
alternative fate that awaits a massive star if its core collapse gen- 
erates a high-energy jet as in HL-GRBs, but in a fashion or within 
the context of a higher-mass stellar envelope that absorbs the jet 
energy before it is able to escape the star. As such, the choked-jet 
events evince no high-luminosity y-ray emission; however, in the 
course of being quenched, their sub-stellar jets may undergo inter- 
nal shocks which could accelerate protons and yield p y and pp 
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neutrinos in the TeV range [46-48]. Observationally, then, the 
choked-jet supernovae would be observed as a relatively low-red- 
shift, highly-energetic supernova (a hypernova) with associated CN 
emission. 

Event rate estimates for choked-jet events are necessarily spec- 
ulative. However, we note that if the process that leads to high-en- 
ergy jet production during core collapse is largely independent of 
the nature of the overlying stellar envelope, then the choked-jet 
supernova rate may be comparable to the rate of HL-GRBs. More- 
over, given that the collapsar mechanism has multiple failure 
modes, including a high-mass stellar envelope, insufficient rota- 
tion, and too much jet precession [49], one can speculate that 
choked jet supernovae should, on generic grounds, be more com- 
mon than HL-GRBs, although potentially (for these and related rea- 
sons) exhibiting less-luminous CN and/or GW emission. 

2.5. Core collapse supernovae 

Core-collapse supernovae (CCSNe) are expected roughly every 
30-100 yr in our Galaxy; given the direct heritage of the revolu- 
tionary neutrino detection of SN1987A [50], any such Galactic cat- 
aclysm would obviously represent a first-grade candidate 
multimessenger transient. 

Seen in relatively nearby galaxies (distances D<100Mpc), 
CCSNe are typically detected by optical observers within a few 
days, although as Swift has demonstrated [51], prompt detection 
of the high-energy shock breakout can be realized with sufficiently 
sensitive instrumentation. On this timescale, the prompt thermal 
(£, > lOMeV) neutrinos should be detectable by IceCube and 
Super-Kamiokande for events within D < 50 kpc [52], 

CN may be produced by these “ordinary” CCSNe, either at shock 
breakout [41 ] or at later times, via interaction of the SN shockwave 
with a dense circumstellar medium [53], The shock breakout CN 
emissions of ordinary CCSNe will be fainter than for the LL-GRBs, 
and so require very nearby events or larger-scale stacking analyses 
to detect with IceCube. CCSNe circumstellar interactions have the 
potential to generate a greater number of CN over an extended 
timescale of weeks to months, and may be detectable with Ice- 
Cube, if candidate CN CCSNe can be identified with confidence 
and in sufficient numbers to make a sensitive search. 

GW emission from generic CCSNe are expected to be relatively 
weak [54], and hence, not detectable except in the case of a Galac- 
tic (or possibly, Small or Large Magellanic Cloud) event. 

Cosmic rays with energies below the “ankle” of the spectrum at 
E os 4 x 10 18 eV are generally attributed to sources within the Gal- 
axy [55]. Since no steady sources have been detected despite sen- 
sitive searches [56], they may be produced in occasional bursts 
within the Galaxy, including Galactic CCSNe. The decay length for 
relativistic neutrons is 9.2 £ 18 kpc, where £ ]8 is the neutron energy 
in EeV. Transient sources of neutrons of 1 EeV or above can thus be 
detected across much of the Galactic disk, including the central 
bulge region, and at higher energies also accessible with Auger, 
sources throughout most of the halo are detectable. 

2.6. Blazars 

Blazars are frequently detected in y-rays at GeV (Fermi, AGILE) 
and TeV (HESS, VERITAS, MAGIC) energies; they are also highly 
luminous X-ray, optical, and radio sources. Since the non-thermal 
X-ray emission of blazars is generally understood to be synchro- 
tron emission from electrons accelerated in shocks of an inner 
jet, shock-accelerated protons are expected as well. The target pho- 
tons for py interactions, leading to CN emission, could then be 
either synchrotron photons from co-accelerated electrons, or opti- 
cal/UV photons originating in the accretion disk or scattered into 
the jet by broad-line region clouds. 


Blazars are highly variable, flaring sources that are much bright- 
er across the EM spectrum during flare episodes than in quies- 
cence. Expected CN fluxes for standard blazar models [21,57,58] 
suggest that typical individual flares cannot be detected with Ice- 
Cube, with the summed contribution of £ 100 such flares required 
to yield an expected >1 CN detection in IceCube. Separately, limits 
on the diffuse CN flux due to the summed contribution of all bla- 
zars in the Northern hemisphere have been published by IceCube 
using the 40-string (roughly one year integration) dataset, con- 
straining some models [24], but without the statistical leverage 
that would be provided by comprehensive EM monitoring yielding 
the times and durations of blazar flares. 

Individual spectral components of blazar flares are also variable, 
and the properties of the highest-energy components are poorly 
constrained due to the difficulty in obtaining sensitive TeV y-ray 
data and simultaneous flare-triggered multiwavelength EM cover- 
age. Moreover, recent Fermi and multiwavelength data suggest 
that spectral breaks and possible additional components are pres- 
ent during some blazar emission episodes [59-61]. Exceptional 
TeV y-ray flares from blazars that lack associated X-ray emission, 
analogs of the 1ES 1959+650 “orphan” TeV flare [62-64], could 
be associated with hadron acceleration that would yield CN fluxes 
well in excess of those expected from typical flares, and hence, 
would be more readily detected by current facilities. 

Energetic GW emission from blazars in the frequency range of 
ground-based detectors is not anticipated. 

2.7. Primordial black holes and other exotica 

If primordial black holes (PBHs) formed in the early universe 
with masses M~5xl0 14 g appropriate for them to undergo 
explosive evaporation at the present epoch [65], then they would 
serve as a distinct and exotic type of multimessenger transient. 

As the black hole loses mass via Hawking radiation, 
dM/dt — -ot(M)/G 2 M 2 [66], its temperature T B h increases, allowing 
an increasing number of particle types (degrees of freedom or dof) 
a(M) to be radiated. Once T B h reaches the quantum chromodynam- 
ic energy scale of Aq C d > 200 MeV, free quarks and gluons will be 
emitted and fragment into hadrons, photons and leptons, resulting 
in a flux of high-energy cosmic ray neutrons and CN [67] that 
might be detectable from Auger and IceCube. 

Detection of PBHs would give dramatic confirmation of Hawk- 
ing’s hypothesis and theories positing early cosmological phase 
transitions (e.g. [68,69]), and enable studies that would likely yield 
deep insights into ultrahigh-energy physics as well as quantum 
gravity. 

Other primordial relics, if they exist, might have decay modes 
yielding harder spectra than the PBHs (e.g. [70-76]). Power-law 
spectra with n oc E 2 or harder would put out a majority of the total 
decay energy at the highest energies, meaning that the first detec- 
tions, or most sensitive limits, on these processes are likely to arise 
from the highest-energy (largest area) facilities, IceCube and Pierre 
Auger. 

3. AMON dataflow and operations 

Multimessenger detection of one of the astrophysical sources 
described above, or of some entirely distinct cosmic phenomenon, 
will require coordination and cooperation between observatories 
of strikingly different design and operation. These observatories 
produce data whose heterogeneity motivates the unified approach 
of AMON for detecting coincidences. Individual facilities participat- 
ing in AMON can be characterized as triggering facilities, follow-up 
facilities, or both. Triggering facilities are typically sensitive to one 
or more messenger type (photons, cosmic rays, neutrinos, or 
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Table 1 

(a) False Positive Rate (FPR) in units of yr 1 for single event streams, including above threshold events that can lead to stand-alone discoveries and the background-dominated 
subthreshold events that are transmitted to AMON. (b) A coarse clustering analysis, shown here for a 'in spatial window and AT = 100 s, will primarily identify pairs of coincident 
subthreshold events with corresponding pairwise FPR. This FPR can be further reduced with a likelihood analysis to tune the distribution rate of AMON Alerts, with the 
subsequent discovery of EM afterglow proving definitive, (c) If three or more subthreshold events are detected by the clustering analysis (shown here for all possible 
combinations that include at least one event from a given stream), or if some other constraint is applied, the FPR is reduced to a level that enables a definitive or near-definitive 
claim of discovery. 




IceCube 

ANTARES 

LIGO-Virgo 

Auger 

BAT 

GBM 

LAT 

HAWC 

(a) Single 

Above thresh. 

-0 

-0 

-0 

-0 

-100 

-250 

-10 

-10 

streams 

Subthreshold 

8.8 x 10 4 
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2.4 x 10 5 
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IceCube 

30 

1.5 

35 
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11 

10 

24 

6.5 


ANTARES 

1.5 

0.5 

12 

1.1 

0.7 

3.5 

7.1 

0.6 


LIGO-Virgo 

35 

12 

N/A 
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53 
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16 

10 

(b) Pair- 

Auger 3 

1.8 

1.1 
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20 
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2.5 

5.9 

1.5 
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BAT 

11 

0.7 

53 

2.9 

N/A 

16 

32 

3.3 


GBM 

10 

3.5 

0.6 

2.5 

16 

N/A 

5.0 
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LAT 

24 

7.1 

16 
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32 
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N/A 

6.8 


HAWC 

6.5 
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10 

1.5 

3.3 

3.2 

6.8 

N/A 


GRB It. curve 5 

0.071 

0.003 

0.16 

- 

0.0004 

0.08 

0.13 

0.019 


SNe It. curve 5 

1.5 

0.07 

3.4 

- 

0.009 

1.6 

2.7 

0.4 

(c) High 

3-fold coinc. 

0.15 

0.03 

0.31 

0.64 

0.12 

0.09 

0.40 

0.08 

significance 

3-fold coinc 3 

0.10 

0.02 

0.15 

0.06 

0.08 

0.04 

0.23 

0.04 


High-sig. EM C 

0.015 

0.002 

0.045 

0.044 

0.010 

0.014 

0.039 

0.005 


PBH search d 

0.13 

0.01 

- 

0.21 

- 

- 

- 

0.35 


a Auger events for the pairwise (and optionally threefold) analysis are selected with galactic latitude within ±5° and energy y 1 EeV. 
b Coincidence rate of event pairs with serendipitous follow-up detection of a GRB or SNe light curve (not including galactic searches). 
c Event pairs with above threshold EM detection. GBM not included as the high significance partner, due to poor spatial localization. 

d An additional temporal cut of AT = 1 s is applied to the pairwise analysis for TeV and higher observatories, to search for primordial black holes and other exotic 
phenomena. Shown here are the FPRs for pairwise coincidences with F1AWC. 


gravitational waves) and search in a wide field of view or monitor 
known sources for transient behavior. From their raw data, they 
construct candidate astrophysical events, here denoted “trigger 
events”, which may transmitted through either public or private 
channels to AMON. These variously represent a single particle 
interaction (e.g., a muon neutrino interaction detected by IceCube), 
multiple detections processed into a high level trigger (e.g., a Swift 
BAT light curve), or a sharp rise in a continuous measurement (e.g., 
an Advanced LIGO event). 

For AMON, a fundamental characterization of each observatory 
is its event-wise false positive rate (FPR). The FPR is the dominant 
component of the total event rate that each observatory sends to 
AMON. The FPR is typically due to intrinsic detector noise or irre- 
ducible backgrounds that produce signal-like events. The FPRs for 
all triggering observatories considered here are listed in Table la. 
Later in this section we describe how these FPRs are calculated. 

By design, AMON will handle data streams dominated by these 
subthreshold signal-like events that individually cannot rise to the 
level of astrophysical discovery. For example, a single muon neu- 
trino detected by IceCube could be of astrophysical origin, but 
the channel is indistinguishable from, and dominated by, atmo- 
spheric muon neutrinos. Similarly, single Galactic neutrons in Au- 
ger cannot be distinguished from the dominant background of 
charged particle cosmic ray events, and lowering the threshold of 
X-ray, y-ray, or gravitational wave detectors increases signal 
acceptance at the cost of simultaneously admitting many addi- 
tional accidental noise events. Flowever, if two or more of these 
events are correlated in arrival direction and time they can attain 
a combined significance that results in AMON issuing an “Alert.” 
An Alert either constitutes an immediate discovery or prompts fol- 
low-up observations that could demonstrate the presence of a 
“smoking gun” electromagnetic counterpart. By having all individ- 
ual observatory FPRs on hand, algorithms running in the AMON 
framework are uniquely capable of calculating the aggregate Alert 
FPR associated with each coincidence detection. The Alert FPR is a 
statistical measure of the quality of a given multimessenger coin- 


cidence, endowing each Alert with a high degree of utility for mul- 
timessenger source searches. 

It is also possible for trigger events from a single observatory to 
be intrinsically above threshold, allowing the observatory to make 
a stand-alone claim for discovery. A burst of three neutrinos in 
close directional and temporal coincidence (IceCube or ANTARES), 
a burst of three air showers in coincidence (Auger), especially if 
their putative source is close enough for them to be neutrons, a sta- 
tistically significant burst of hard X-rays (Swift BAT, Fermi GBM) or 
y-rays (Fermi LAT or HAWC), or a strong strain registered by a 
gravitational wave detector (LIGO-Virgo) would each constitute a 
discovery. For these stand-alone discoveries, AMON can serve as 
a quick and convenient conduit for disseminating the source coor- 
dinates, enabling timely follow-up observations at multiple dispa- 
rate observatories. AMON will leverage the GRB Coordinates 
Network (GCN) [77], among others, to promulgate its Alerts. 

Underpinning AMON’s sensitivity to new phenomena is the 
broad (multiple steradian) sky coverage and high duty cycles of 
the triggering observatories. Fig. 1 shows the results of a simula- 
tion of a calendar year for the triggering facilities that have signed 
the AMON MOU (IceCube, ANTARES, Auger, HAWC, and Swift BAT) 
and potential future signatories (Fermi LAT and GBM, and LIGO- 
Virgo), establishing that - absent severe disruptions - at least 
two facilities are always observing any given part of the sky, and 
that about 94% of 47t sr-yr lies within the field of view of three or 
more facilities. 

Existing multimessenger searches involving only pairs of obser- 
vatories are typically in a master-slave relationship, e.g., a trigger- 
ing observatoiy initiates follow-up at an optical telescope [78-80]. 
In contrast, AMON will enable archival and real-time searches for 
coincidences among multiple observatories in a peer-to-peer 
relationship that will provide markedly increased opportunity for 
discovery of new multimessenger sources, as detailed in Section 4 
below. The real-time AMON analysis will be informed by the 
extensive effort already invested by the triggering observatories 
[78,81-90]. The trigger events these observatories produce will 
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Fig. 1 . (left) Distribution of the total number of trigger facilities observing a cosmological source, averaged over time and sky location, (right) Average number of 
observatories simultaneously viewing a source, as a function of location. Spatial variance is caused by the limited field-of-view of some observatories and the movement of 
orbital telescopes. 


be heterogeneous, covering the full range of messenger particles 
and differing significantly in background rate and localization pre- 
cision. In some cases the events will include additional quality 
information that may be useful in refining the single-stream false 
positive rate. For purposes of AMON, all events will be couched 
in the common language of statistics, containing information 
including a trigger time, event position, and a positional error or 
probability density function (PDF). 

Once the individual event streams have been characterized, 
including tracking observatory pointing and FPR, the multiple 
event streams can be combined. At the most basic level, this will 
be carried out via a coarse clustering analysis, searching for tempo- 
ral and spatial coincidences. The output of the clustering analysis is 
dominated by pairwise coincidences, with the FPR for an ordered 
pair of triggers from observatories (a, b) given by, 2 

€ P) ~*t a k AT{n * b)AQab ’ (1) 

where R a is the subthreshold event rate from observatory a to 
AMON (assumed constant here), Q a is the field of view of observa- 
tory a, AT is the total width of the temporal search window, ( Cl ab ) is 
the time-averaged overlapping field of view of observatories a and 
b, and A Cl ab is the search area for coincidences between observato- 
ries a and b. In Table 1 (row “b”), we show the result of this calcu- 
lation for AT = 100 s and Afi a /, adjusted to give ~90% acceptance for 
signal event pairs (where each event is at the threshold for trans- 
mission to AMON). The self-coincidence of the individual observa- 
tories is given by the diagonal R aa (although this calculation is not 
appropriate in the case of preprocessed event streams) and the 
sum R ab + R ba (a ^ b) gives the total number of pairs from a and b 
without respect to order. As expected, the combined FPR due to a 
pair of false positive events is orders of magnitude lower than that 
of the individual component subthreshold event streams. 

This particular calculation is by no means definitive. A trivial 
adjustment of the cuts will increase the rates (useful for generating 
AMON Alerts for testing purposes). More importantly, additional 
information can drastically lower the FPR for the AMON analysis 
(Table 1, row “c”). For example, one can statistically stack the 


2 In some cases, it may be advantageous to use one event stream to trigger event 
selection for another, e.g., if the rate of single events is high, as with HAWC, or if the 
localization is poor, as with low-significance gravitational wave events. However, 
processing of an observatory’s data can always be carried out, either at AMON or prior 
to transmission, to yield a manageable single-observatory FPR, and such preprocess- 
ing has been assumed for this discussion. 


Alerts though an archival study, or shorten the temporal search 
window to search for exotic phenomena. In real-time, one can re- 
quire that at least one event be of high significance or that the 
coincidence contains at least three events. 3 

A promising feature of AMON is its ability to distribute Alerts in 
real-time to initiate follow-up observations. The relatively high 
FPR represented by the sum total of all pairwise coincidences 
among the observatories considered (Table 1, row “b”) may be tol- 
erable for some follow-up facilities, especially if their response is 
highly automated. However, many follow-up facilities will have 
limited observing time available for this science and may require 
a substantially lower FPR for triggering. This can be achieved via 
a refined analysis, providing a combined best fit location and error, 
as well as a maximum likelihood ratio or Bayesian probability. The 
likelihood or probability measure is drawn from a continuous dis- 
tribution, in contrast to the FPR of the clustering analysis which, gi- 
ven a fixed event rate for each stream and a fixed set of cuts, takes 
on a unique value for any pair of event streams. Prior information 
can easily be incorporated into the likelihood analysis including, 
for example, previous detection limits, observatory sensitivity, or 
a galaxy catalog, and these choices can be tuned for different fol- 
low-up programs. 

As an independent channel of scientifically-relevant data and 
constraints, follow-up observations have the potential to provide 
the added information required for a definitive discovery. By distrib- 
uting candidate source positions in real-time to follow-up facilities, 
with statistically valid measures of the FPRs, AMON will enable effi- 
cient fast-response counterpart searches and studies across the elec- 
tromagnetic spectrum. The existence and nature of such EM 
counterparts may prove decisive in verifying the existence of some 
of the first multimessenger sources, and can reveal the nature and 
detailed properties of the source, including important ancillary 
information such as the source environment, host galaxy (if any), 
and properties of intervening gas and dust along the line of sight. 
The discovery potential of pairwise coincident events with positive 
follow-up detection is described in more detail below. 

AMON will be deployed in three phases. During the first year of 


3 The FPR for a coincidence of the ordered triple (a. b, c) is given by 
^abJ ~ 2 o)o(a (ATAfl aftc ) 2 (Qafljc), where the factors AQ abc , an equivalent search area 
for a threefold coincidence, and (fi of , c ), the joint field of view, are determined by 
Monte Carlo. The factor of j depends on how the temporal search is defined, where 
here we require that all three events occur within AT. In Table lc, we present the 
totals ■ for the case where all Auger data is included and for the case where 

Auger data is restricted to the galactic plane. 
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operation beginning in mid-2013, AMON will use archived data to 
develop and tune analysis algorithms. Such archival analyses, while 
a necessary first step prior to activation of real-time alerts, will also 
be of intrinsic scientific interest, and may point to possible new 
astrophysical signals or new and constraining upper limits on 
jointly-emitting source populations. Even as AMON moves toward 
real-time operation, certain source populations will be best discov- 
ered through a collaborative blind study of accumulated data; for 
example, in searches for exotic particles, as discussed below for pri- 
mordial black holes. To enable this, participating AMON members 
will, at any time, be able to utilize the AMON database and data prod- 
ucts, under rules that are established in the AMON MOU [91]. 

In subsequent operational phases, standard AMON analyses will 
be modified to operate in real-time, enabling identification of can- 
didate transient sources as soon as the events are transmitted. At 
first, the effort will focus on those discoveries that can be made 
by the triggering facilities alone, with FPRs similar to those in 
Table lc. Even at this stage, the real-time nature of the network 
plays a critical role, providing the triggering facilities with AMON 
Alerts that can prompt a closer analysis of their own data. The 
Alerts shared at this stage will include many whose combined 
FPR is too high to claim detection of a transient. These will not only 
exercise the AMON system, but may initiate a deeper search by 
those participating observatories that did not initially report a trig- 
gering event. An example of this discovery mode is the joint search 
for high-energy cosmic neutrinos and gravitational waves, as dis- 
cussed below. 

As quickly as is feasible, AMON will transition to its final phase, 
with Alerts distributed to participating follow-up observatories to 
enable the near real-time search for electromagnetic counterparts. 
For sources that exhibit detectable afterglow, this is the most pow- 
erful technique, leveraging the high acceptance of the initial FPR 
threshold, yet resulting in a final FPR that is lower than that avail- 
able by the other techniques. The estimated final FPRs are shown in 
Table lc, for the case of serendipitous follow-up detection of GRB 
afterglow (assuming Swift XRT sensitivity and a ~1° search region) 
or a SNe light curve (assuming the detection limits of [78]). In the 
next section, as a proof of concept, we show that significant sensi- 
tivity gains are made by AMON in the search for electromagnetic 
counterparts to cosmic neutrino sources. 


4. Discovery power of AMON 

In this section we attempt to quantify the discovery power of 
AMON. Any such effort must grapple in some way with the vast ar- 
ray of possible discoveries that AMON might contribute to, both in 
terms of the various candidate source populations (Section 2) and 
in terms of the diverse set of possible partner facilities (Section 3), 
including the combinatorial implications of joining those facilities 
together into a single network. 

We take two distinct approaches here. In the first approach, we 
focus on three specific examples of multimessenger source popula- 
tions, construct models for the populations and their multimessen- 
ger signals, and simulate the performance of appropriate partner 
facilities both with and without the contributing capabilities of 
AMON. Specifically, Section 4.1 examines how the distribution of 
AMON Alerts can greatly enhance the search for electromagnetic 
counterparts to candidate sources of cosmic neutrinos, while Sec- 
tion 4.2 quantifies the gains for a search for joint sources of gravi- 
tational waves and CN. Section 4.3 discusses an example search for 
exotic phenomena, showing how an AMON analysis may provide 
the first clear signature for primordial black holes. 

In the second approach, we explore the sensitivity of the various 
possible AMON partner facilities in a relatively model-independent 
fashion (Section 4.4). The diverse array of messenger types, above- 


threshold and sub-threshold trigger definitions, and likely source 
populations (e.g., Galactic, Local Group, Extragalactic, cosmologi- 
cal) addressed simultaneously in this analysis make an assump- 
tion-free exploration impossible; however, we have adopted a 
relatively simple set of default assumptions that at least enable 
the cross-facility comparison to be made. In doing so, we address 
the merits and demerits of this approach as compared to our pre- 
ferred simulations-based method. 

4.1. Follow-up of candidate cosmic neutrino sources 

As a first example, we consider the follow-up imaging of candi- 
date cosmic neutrino sources. Such searches are already underway 
via efforts by individual neutrino observatories and their follow-up 
partners. For example, in [92], the IceCube Observatory identifies 
clusters of two or more neutrino events as a trigger for follow-up 
imaging. Due to the high rate of atmospheric neutrinos and other 
backgrounds, the expected FPR for a pair of up-going events for Ice- 
Cube-86 is approximately 30 yr -1 (Table 1, row “b”). Alerts are 
then distributed to the ROTSE telescope array or, after using a max- 
imum likelihood analysis to refine their number, to Swift and PTF 
[80,93], This program is already producing the first limits on SNe 
with choked jets [78]. 

As a shorthand, we introduce the notation “2v” to indicate a 
trigger condition of 2 or more coincident neutrinos, or v-Ny for 
one or more neutrinos and N or more y-rays. To be concrete, we 
study the follow-up imaging of 2v alerts at a rate of 10 yr -1 using 
Swift’s co-aligned X-ray and UV/Optical Telescopes (XRT, UVOT), 
with this considered a status quo approach for the scientific com- 
munity. Due to the narrow field of view (0.4" diameter) of the 
XRT and UVOT, 7 Swift pointings are required to cover the error 
region of each alert, totaling 70 pointings in all. 

In contrast, the AMON approach enables the realtime coincident 
analysis of neutrino and electromagnetic data streams, comparing 
the arrival time and direction of a single candidate CN with data 
from Swift’s large field of view Burst Alert Telescope (BAT), Fermi’s 
Large Area Telescope (LAT), or the upcoming ground-based HAWC. 
Because of the superior AMON localizations, derived via the joint 
analysis with high-energy y-rays, the same follow-up telescope re- 
source (70 pointings) can be used for follow-up of a significantly- 
increased number of AMON alerts. Here we consider a program 
that follows up five 2v alerts (at 7 pointings each) and 35 v-Ny 
alerts (1 pointing each). For the v-Ny alerts, we considered the 
pairing of IceCube with both BAT and LAT. 

The relative sensitivity of the status quo and AMON approaches 
was studied via a Monte Carlo calculation. As a source model, we 
assumed that both y-ray and neutrino spectra follow a broken 
power-law with initial and final logarithmic slopes of a = 1 .0 and 
P = 2.0, and break energies g y = 0.2 MeV and = 0.35 PeV, 
respectively. The neutrino and y-ray fluence was allowed to vary, 
as indicated by the axes of Fig. 2. Using realistic models for the Ice- 
Cube-86, BAT, and LAT effective areas, point spread functions, and 
their overlapping field of views, we were able to determine the 
acceptance of each simulated source. Realistic background models 
were adopted for each observatory, enabling the calculation of 
false positive rates for each joint analysis. 

In each case, the FPR was tuned via a likelihood analysis of the 
relative position of the v and y events. Subthreshold BAT and LAT 
events were allowed into the simulated data stream, with terms 
added to the likelihood function to favor the more significant EM 
signals. Furthermore, we added a regulating term to the likelihood 
function to account for the possibility of a real y-ray source in acci- 
dental coincidence with an IceCube background event (e.g., an 
atmospheric neutrino). Once the FPR is chosen for any pair of 
observatories, we identify the required likelihood threshold, and 
can then calculate the corresponding sensitivity to our modeled 
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Fig. 2. Sensitivity of the AMON multi-channel method compared to the status quo, 
presented as a ratio ranging up to ~10 6 , for follow-up observations of candidate 
astrophysical neutrino sources. The assumed shapes for both gamma and neutrino 
spectra follow a broken power law with initial and final logarithmic slopes of 
at =1.0 and (1 = 2.0, and break energies g 1 = 0.2 MeV and 7;\. = 0.35 PeV. The 
vertical line shows a typical threshold for a hard X-ray observatory (e.g., Swift BAT), 
and diagonal lines corresponding to predicted neutrino to gamma-ray fluxes under 
a range of current theoretical GRB models [22], 


source. While the initial FPR is quite high (10-70 alerts yr _1 ), the 
final FPR will be orders of magnitude lower after the inclusion of 
follow-up data and, as such, can be ignored in the study of sensitiv- 
ity. While here we have characterized the sensitivity of the trigger- 
ing methods (AMON versus status quo), there are a number of 
unknown factors in the sensitivity of the subsequent follow-up 
searches; for example, due to optical magnitude. However, we 
study the ratio of the sensitivity of the two methods, expecting 
the acceptance of the two follow-up programs will be approxi- 
mately the same and to cancel in the ratio. 

The results are shown in Fig. 2, where substantial gains in sen- 
sitivity (up to 10 6 ) are realized. The sensitivity of each method is 
defined as the fraction of true sources pursued using a given fol- 
low-up capability. In addition to enabling pursuit of an increased 
number of alerts overall, the AMON approach realizes these large 
sensitivity gains by triggering from single neutrinos rather than a 
neutrino pair, leveraging the increased interaction strength of y's 
over v’s. The left axis of Fig. 2 is defined in such a way that it could 
equally well describe the expected event rate for ANTARES, with 
only a small change in the contours due to the narrower neutrino 
point spread function. 

Gains of this type have also been demonstrated for gravitational 
wave source searches, using GW-y triggering over the status quo 
GW-only methods [94]. For GW-only searches, the poor GW locali- 
zation is filtered by requiring coincidence between the GW localiza- 
tion and a nearby cataloged galaxy [79,95], While useful, this 
approach suffers from several shortcomings: First, it depends on 
the completeness and astrophysical relevance of the associated gal- 
axy catalog; second, it cannot be extended to the subthreshold re- 
gime due to the dramatic increase in false alarm rates; and finally, 
the approach is limited to an horizon of <100 Mpc, beyond which 
galaxy fields become crowded and catalogs are highly incomplete. 
The multi-channel, multimessenger approach by contrast provides 
a unique candidate source position for follow-up pointing, and lever- 
ages the coincidence technique to push into the subthreshold, possi- 
bly revealing otherwise-undetectable source populations. With 
Advanced LIGO expected to ultimately achieve sensitivity out to 
> 500 Mpc [96], the AMON approach can offer a >125-fold gain in 
sensitivity for events with bright EM counterparts, or >1 000-fold 
for regions that have the least complete galaxy catalog. 


Table 2 

Effective volumes for a joint GW-CN search, assuming = 10 2 M,. and 
E“° = 10 49 erg for a choked GRB model. 


Analysis 

Integration limit 

Neutrino 

yitff) 


(Mpc) 

threshold 

[Mpc 3 ] 

GW+CN 

80 

n v >3 

4.6 x 10 2 

AMON 

160 

n v ^ 1 

5.9 x 10 4 


4.2. Joint detection of neutrinos and gravitational waves 


Here we consider the joint detection of gravitational waves and 
high energy cosmic neutrinos as an example of a science investiga- 
tion enabled by the realtime operation of AMON. The search for 
jointly emitting GW+CN sources has been addressed in [97,98], 
including predictions for Advanced LIGO-Virgo and IceCube-86, 
and we follow the same general approach (see also [99,100]). How- 
ever, while previous authors focus on the joint analysis of a multi- 
messenger data archive, we emphasize the possibility of carrying 
out this analysis in realtime. We characterize the sensitivity of a 
blind (untriggered) search by Advanced LIGO-Virgo using the hori- 
zon distance D c w. inside which a GW source is expect to generate 
an above-threshold detection, 


Dew 63 80 Mpc 


( 2 ) 


where is the isotropic equivalent gravitational wave energy of 
the source. This is based on the results of [101] and assumes little 
of the GW waveforms. Without AMON-like infrastructure, a real- 
time detection of both GW and CN requires that a source be within 
this horizon and that it produce a neutrino flux sufficient to trigger 
n v ^ 3 events in IceCube or ANTARES (since n v ^ 2 would lead to an 
inconclusive false positive rate). These detections would then be 
compared post facto, either through an archival analysis or above- 
threshold alert distribution (e.g., through GCN), to determine if there 
are any joint GW-CN sources. 4 

By contrast, a system like AMON will provide two-way commu- 
nication between observatories, enabling a realtime triggered 
search of GW data, either by transmission of pre-alerts from AMON 
to LIGO-Virgo or by carrying out a triggered search on subthreshold 
GW data transmitted to AMON (albeit in some preprocessed form). 
Searches of this kind suffer a lower trials penalty and may assume 
various templates for the sources. For example, the postulate that 
CN are detected in coincidence allows one to make an assumption 
about the orientation of the source, making the GW signal more 
favorable for detection [102]. In all, we estimate that these effects 
will combine to increase the horizon distance to 2 x D c w [102]. 
Furthermore, only a single neutrino event is required to trigger 
the analysis, providing significant gain for this approach since the 
expected number of neutrinos has been observed to be typically 
« 1 [103], 

Since we are considering the relatively deep horizon distance 
for Advance LIGO-Virgo (so that galaxy distribution can be as- 
sumed uniform), but not so far that redshift effects become impor- 
tant, it is sensible to apply the effective volume approach of A, 
integrating up to the physical limit of D c w or 2D C w- The resulting 
effective volumes are given in Table 2. 

Under the assumption of a low false positive rate, the 90% upper 
limit (chosen for consistency with [97]) for the joint GW-CN source 
density is given by 


4 In [97], the authors consider a slightly different problem, being interested in the 
above-threshold detection of GW or CN. While an interesting measure, their 
calculations do not provide the sensitivity for above-threshold detection of joint 
sources that we need for comparison to AMON. 
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Fig. 3. (Left) Projected GW+CN source population upper limits after one-year of AMON utilizing observations with Advanced LIGO-Virgo and IceCube-86. (Right) Ratio of the 
sensitivity for the AMON approach to the sensitivity of the status quoa pproach (as described in the text). 


Pgw-cn < 


2-3 f„ 

jy(eff) 


[yr- 1 Mpc- 3 ], 


( 3 ) 


where T is the livetime, which we take to be one year, V (eff) is the 
solid angle averaged effective volume for the given method, and f b 
is the beaming factor for neutrinos. We take the latter to be 
f b — 14 (from [97]) but it will cancel in our final result. 

The effective volume depends on through Eq. 2. For neutri- 
nos, the expected number of detections from a fiducial source at 
distance r can be related to the neutrino isotropic energy E 1 ™, using 
a linear relationship provided in [97], 


n v 


\( E' so \/ 10Mpc \ 2 
K (lO 49 ergy [ r J ’ 


( 4 ) 


where, for IceCube-86, K = 0.75 for a high luminosity GRB model 
and k = 1.6 for a model of choked GRBs. 

The results are shown in Fig. 3. The left side shows the upper 
limit in units of yr- 1 Mpc~ 3 that would be obtained by a non-detec- 
tion with AMON. The right side shows the multiplicative gain in 
sensitivity that the AMON approach provides over an analysis that 
requires above-threshold detection in both GW and CN channels. 
The gain is most significant for bright GW sources and weak neu- 
trino sources, where the burden for discovery is shifted from the 
CN to GW channel. However, the gain is strictly greater than 1, 
with an asymptote of 8, representing the increase in effective vol- 
ume due to triggering the GW search. 


4.3. Searching for primordial black holes 

Here we consider the search for primordial black holes as an 
example of a collaborative study of accumulated multi-facility 
data, as enabled by the AMON infrastructure and data sharing 
policies. In particular, we have modeled a joint search with the Ice- 
Cube neutrino observatory, HAWC TeV gamma-ray observatory, 
and the Pierre Auger cosmic ray observatory. All are sensitive to 
particles produced in the final stages of PBH evaporation, including 
a potential signal from ultrahigh-energy neutrons detected at 
Pierre Auger, since these neutrons would not suffer from the mag- 
netic deflection and time delay effects of charged cosmic rays 
(the chief UHE neutron background). 

Based on the energy threshold of these experiments and energy 
dependence of various models predicting energy spectra of PBHs in 
their final stage of evaporations (e.g. [104-106]), the model of a 
non-rotating, uncharged black hole without a chromosphere by 


MacGibbon and Webber [104] was chosen for estimation of the ex- 
pected signal at each of these detectors. The time integrated parti- 
cle spectra above 100 GeV were computed following methods from 
Refs. [107,104,108,109], The main particle decay chains considered 
in our calculations are described in B. 

It is anticipated that, when fully operational, HAWC by itself 
will either detect or provide the best upper limit on direct searches 
for PBH decay, around 2 orders of magnitude better than the cur- 
rent upper limit (as shown in Fig. 4). To do so, HAWC must perform 
a blind search of its own data. After a cut to remove many hadronic 
events, 5.1 kHz of subthreshold trigger events remain, which will 
be analyzed for coincidences across a spatial template with 10 4 tri- 
als [110]. Taking a generous time window of AT = 5 s and a num- 
ber threshold of n 3 : 20, the HAWC false coincidence rate can be 
reduced to <1 yr~\ However, it is expected that HAWC will also 
observe real astrophysical transients that are unrelated to PBHs, 
forming a second type of false positive for the PBH search. We con- 
servatively estimate this number to be R[ FP) <17 bursts per year, 
approximately equal to the known LAT GRB rate, allowed to fluctu- 
ate up by 2a. If, after T = 1 yr of livetime, 17 unidentified bursts 
were observed, the 99.73% (3er) upper limit for PBH would be given 

by 


32.1 

Ppbh < 


[yr 1 pc 3 ] , 


( 5 ) 


where 32.1 is chosen to give the 99.73% upper limit of the Poisson 
distribution, and V (eff) is the integral of the trigger condition 
P(n > 20|r) over 4nr 2 dr, averaged over solid angle, where r is the 
distance from the Earth. The value of V (efF) is specific to the sensitiv- 
ity of the observatory and the intrinsic y-ray fluence of the fiducial 
source, and is further described in A. For the current calculation, we 
estimate V (0ff) = 0.037 pc~ 3 . The calculated value of p PBH <3.5x 
10 3 yr~'pc~ 3 can be compared to the current best limit of 
1.9 x 10 5 yr- 1 pc- 3 from VERITAS [111], 

Between the current VERITAS limit and the projected HAWC 
limit, there exists the possibility of a future positive detection of 
PBH decays. However, the HAWC data by itself may be unable to 
distinguish between PBH and other scenarios. A joint AMON study 
of HAWC data with neutrino and neutron data (modeled here with 
IceCube and Auger), will enable the search for multimessenger 
coincidences, providing a unique PBH signature. Just a few such 
coincidences should reveal a timing and energy structure that are 
indicative of Hawking evaporation. 
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Fig. 4. Expected number of sources detected in at least two of the v,y, and n 
channels, modeled here with the lceCube-86, projected HAWC, and Auger sensi- 
tivity, as a function of the primordial black hole density. The shaded region to the 
right indicates PBH densities that previous direct searches have excluded, while the 
vertical line shows where the upper limit is expected to move if HAWC does not 
observe a gamma ray excess after a year of livetime. In between, there is the 
possibility of positive joint detections, with potentially a unique PBH signature. 


Table 3 

Effective volumes for multimessenger PBH search. Values are averaged over 47t sr, 
with only those regions of overlapping sensitivity contributing to each pairwise 
calculation. 


Observatories 

Trigger conditions 

Ipc 3 ] 

vSS- IPC 3 ] 

A 

B 

A 

B 



HAWC 

- 

n y ^ 20 

- 

0.0374 

0.245 

IceCube 

HAWC 

n v ^ 1 

tly ^ 13 

9.8 x 10- 5 

6.4 x 10~ 4 

IceCube 

HAWC 

n v ^ 2 

n 7 > 13 

1.3 x 10~ 6 

8.3 x 10~ 6 

HAWC 

Auger 

n y > 13 

n„ » 1 

1.3 x 10~ 7 

8.6 x 10~ 7 

Auger 

IceCube 

n n ^ 1 

n v ^ 1 

4.9 x 10~ 9 

3.2 x 10~ 8 


Folding together the PBH source model from B and the volume 
integrals of A, we calculate the effective volumes given in Table 3. 
To maximize sensitivity to neutrinos and neutrons, we use number 
thresholds of n v ^ 1 and n„ > 1. For y-rays, the low false alarm 
rate can be relaxed significantly by requiring coincidence with a 
single v or n, allowing the HAWC threshold to be lowered to 
n y Ss 13. The expected rate of true positive joint detections in chan- 
nels a and b is then given by 


4.4. Cross-facility sensitivity comparison 

In general, our preferred approach to evaluating the discovery 
power of AMON has been via direct simulation of modeled source 
populations as in the sections above. The chief shortcoming of this 
approach is that, by focusing on individual case studies, it fails to 
provide a comprehensive overview of the discovery space that will 
be opened up by the array of multimessenger facilities now coming 
online - and by AMON itself. 

For completeness, then, in this section we present a relatively 
model-independent cross-facility comparison of the sensitivity of 
candidate AMON partner facilities. 

Before presenting this analysis, however, we wish to point out 
some key limitations. First, quantifying the discovery power of a 
multimessenger search is intrinsically and strongly model-depen- 
dent. For example, the relative probability of triggering a high-en- 
ergy y-ray observatory compared to a low-energy X-ray 
observatory depends strongly on the spectral shape of the source, 
with a steeply-falling spectrum favoring the low-energy facility, 
and relatively hard spectra favoring the high-energy facility. 

Second, the discovery power of AMON depends critically on the 
joint multimessenger emissions of entire populations of sources. 
For example, jointly emitting vy sources (e.g., HL-GRBs) may be 
routinely detected in gamma-rays and undetected in neutrinos. 
However, given a large population of such sources (e.g., hundreds 
of HL-GRBs per year with current satellites), each with low neu- 
trino fluences, occasionally one source will yield detection of a sin- 
gle neutrino. Modeling the source population as in Section 4.1 
allows us to take this “rare but significant” scenario into account 
when calculating the increased discovery power of AMON for 
jointly-emitting vy sources. 

Finally, the ability to make joint detections via two or more 
messenger types depends critically on the relative energy fluences 
in each channel. Yet the fluences in different channels are not gen- 
erally expected to be equal, and in some cases are not even ex- 
pected to be well-correlated (except for the correlation induced 
by observing sources across a substantial range in distance). In 
general such effects can only be addressed realistically via 
source-specific simulations. For example, the total energy emitted 
in very high energy (VHE) neutrinos and y-rays by HL-GRBs (and 
potentially LL-GRBs) is expected to be comparable in some models; 
however, the VHE y-ray emissions will be strongly attenuated by 
interactions with the extragalactic background light beyond a rel- 
atively short cosmic horizon of D ~ 1 Gpc. 

For purposes of our cross-facility comparison, we ignore these 
effects and choose a single standard source spectrum, a power- 
law with spectral slope dN/dE oc E~ 2 . This is the same spectrum 
as used in Section 4.1, but different from the spectra used in our 
other examples. The spectrum is extrapolated without any break 
to lower energies when needed, which also differs from our ap- 
proach in the source specific examples. The chief advantage is that 
it enables comparison of observatories operating at radically differ- 
ent energies using the single parameter F y , defined via 


R 


(TP) 

ab 


Ppbh 

C’ 


[yr 1 pc“ 


(6) dN y _F y E- 2 
dEdA ~ TnTO ' 


( 7 ) 


which we show as a function of p PBH in Fig. 4. Over the range of p 
which might lead to future positive detections by HAWC, the total 
number of expected multimessenger coincidences ranges from 0.1 
to 100 yr _1 . Most of these are v-Ny (i.e. a coincidence between a sin- 
gle neutrino and N y-rays). However, there is a non-negligible pos- 
sibility for observing n-Ny coincidences (where the neutron is 
detected by Auger) or 2 v-Ny, providing an exciting possibility for 
detecting three subthreshold events in coincidence (where the Ny 
from HAWC are treated as a single subthreshold event). 


Defined in this way, F-. is the gamma energy fluence (eV cm~ 2 ) per 
decade of energy. 

For neutrinos and neutrons, we use the same E~ 2 shape and re- 
place the normalization parameter F y with F r or F„, as appropriate, 
for the neutrino or neutron energy fluence per decade. We note 
again that EeV neutrons will only be detectable for sources within 
the Galaxy. The various parameters F r ,F v , and F n can then be rep- 
resented on a single axis in the same figure. However, naturally 
these remain entirely different quantities, representing the energy 
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Fig. 5. Conditional triggering probability P as a function of energy fluence F (with 
linear scale for P 0.5 and log scale for P < 0.5), assuming the source is in the field 
of view and the observatory is active. For BAT, LAT, HAWC, IceCube, and Auger, the 
bottom axis shows the fluence per decade under the assumption of an F 3 spectrum 
while, for the Advanced LIGO model, total fluence is shown. 

fluence per decade at Earth for the associated messenger only, 
without any implication that likely source models are expected 
to emit the different messenger particles in equal or fixed propor- 
tion. Completing a characterization of cross-facility sensitivities 
would require adjusting for the emission properties of a specific 
source population, as in our simulations, and accounting for any 
propagation effects that would affect the fluences observed at 
Earth. 

Finally, we can incorporate gravitational wave detection into 
the same scheme by following the equations of Section 4.2 and 
using the shape of the efficiency curves from [101]. We assume 
that a total GW energy of O.OlM 0 c 2 will produce a GW event at 
the nominal detection threshold when located at 160 Mpc from 
Earth (assuming coincidence with some other messenger type in 
AMON). We note that the total event fluence F C w derived in this 
way is fundamentally different from the F y ,F v , and F„ fluences, 
which are defined as fluences per decade. At the same time, we note 
a strong expectation that, for the strongest GW transients, the GW 
energy release will be orders of magnitude greater than for any 
other messenger type. 

The result of these calculations is presented in Fig. 5 as a plot of 
the conditional trigger probability P for each AMON partner facility 
versus the event fluence or fluence per decade for the appropriate 
associated messenger, as defined above. The conditional trigger 
probability is the probability of detection for an event occurring 
within the facility field of view while the facility is in active oper- 
ations at nominal expected sensitivity. 

In order to derive a single conditional probability of triggering 
for any event of particular fluence, we have made some further 
assumptions. With respect to the y-ray facilities BAT, LAT, and 
HAWC, we assume that the energy fluence arrives in one second 
and that the observatory is set to trigger on one-second bursts; 
as a trigger threshold, we require that the expected number of 
event photons exceed the expected background in this interval 
over a region sized to the facility PSF. This approach leads to a stee- 
ply-rising trigger probability as the threshold condition is ap- 
proached and exceeded. Moreover, we note that the ordering of 


sensitivities shown in the plot for BAT, LAT and HAWC depends 
strongly on the spectral model. Had we chosen to break the spec- 
trum, making it flatter at low-energy, then LAT and HAWC would 
have been more competitive relative to BAT. 

For each observatory, we can identify a nominal threshold cor- 
responding to P = 50% conditional triggering probability, which 
we see occurs at F y < 10 5 eV cnr 2 for BAT, LAT, HAWC. In contrast, 
IceCube has a greater nominal threshold of F„ ~ 10 8 eV cnr 2 . 
Among events with comparable y-ray and neutrino fluences (at 
Earth), only the very brightest y-ray sources would be expected 
to yield enough neutrinos to exceed this threshold. However, as 
discussed above, an occasional neutrino detection can also be 
anticipated as the result of statistical fluctuations among a much 
larger population of weaker sources. By virtue of AMON coinci- 
dence analysis, this neutrino detection can be assumed to be effec- 
tively background-free, which results in a sensitivity curve that 
rises much less steeply than that of the y-ray observatories. Specif- 
ically, at F v = 10 5 eV cnr 2 , two or three of 10 3 sources would be 
detectable in neutrinos, a possibility that is highly relevant both 
for GRBs and for possible unknown source populations. 

With respect to the Pierre Auger Cosmic Ray Observatory, Fig. 5 
apparently shows a “best of both worlds” situation, with a rela- 
tively low nominal threshold of F„ < 10 5 eV cnr 2 and a gently fall- 
ing trigger probability below that. This has the same shape as 
IceCube’s P(F) curve, given by the Poisson probability of contribut- 
ing at least one event to a coincidence analysis, but with a higher 
event yield for a given fluence. However, it should be recalled that 
the Auger calculation is for neutron events only, which probe a 
highly limited detection volume due to their half-life. Furthermore, 
the assumed FT 2 spectrum cannot extend indefinitely to high en- 
ergy lest it require infinite energy at the source; steeper spectra, 
or spectra with breaks within or below the Auger range, would re- 
duce the relative sensitivity of Auger by comparison to the other 
facilities. 

We conclude that, while a unified presentation of cross-facility 
sensitivities can be achieved (Fig. 5), the results retain a significant 
number of assumptions such that interpreting the figure, and 
translating it into real-world expectations for actual cosmic events, 
is not straightforward. 


5. Summary and conclusions 

The Astrophysical Multimessenger Observatory Network under 
development at Penn State will link multiple high-energy, multi- 
messenger, and follow-up observatories together into a single 
comprehensive system. We have explored the scientific opportu- 
nity for AMON (Section 2), which is centered on the discovery 
and exploitation of multimessenger transients. The nature of the 
brightest such transients, which may manifest as the first detec- 
tions of gravitational waves, high-energy cosmic neutrinos, or 
high-energy cosmic ray neutrons, is still unknown; likely candi- 
dates include blazar flares and a diverse array of y - ray bursts 
and supernovae, as well as exotic phenomena such as the evapora- 
tion of primordial black holes from the early universe. 

We have described the design, infrastructure, and current and 
projected partner facilities of AMON (Section 3), showing that 
the wide fields of view, high duty cycles, and subthreshold event 
rates of the facilities are such that a robust and automated statis- 
tical search for coincident events seen in the data of two or more 
facilities is both interesting and feasible (Table 1). In particular, 
false positive rates for the resulting AMON Alerts are sufficiently 
low that comprehensive ground-based optical follow-up cam- 
paigns can be contemplated. 

With neutrino and y-ray observatories are already participating 
in AMON, and discussions underway with cosmic ray and 
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gravitational wave observatories, as well as multiple follow-up 
observatories, AMON is poised to begin real-time operations with- 
in a year. 

To demonstrate the power of the AMON approach, we have car- 
ried out three sets of simulations using current theoretical models 
of multimessenger phenomena, and compared the sensitivities of 
the various multimessenger facilities to one another using a some- 
what simplistic ansatz (§4). In the first simulation, we explored the 
challenge of EM counterpart searches for candidate GRB-associated 
CN observed in IceCube, showing that vetting candidate CN events 
against multiple EM data streams realizes a >1 000-fold gain in the 
efficiency of EM follow-up observations (Fig. 2). In the second sim- 
ulation, we explored the improvement in search sensitivity for 
jointly emitting GW+CN transients that is realized by extending 
these searches into the subthreshold regime, rather than restrict- 
ing the search to events that generate statistically-significant sig- 
nals in both channels; we find that a >1 0-fold increase in event 
rates, or improvement in upper limits, is easily achieved (Fig. 3). 
In the third simulation, we explored the multimessenger signature 
of PBH evaporation, demonstrating that - consistent with current 
upper limits on the local PBH density - the coincidence of an Ice- 
Cube-detected neutrino with a cluster of HAWC-detected y-rays 
could provide evidence for observation of a PBH evaporation event 
within the first year of HAWC operations (Fig. 4). The promise and 
challenges of multimessenger transient detection seen in these 
three case studies are also illustrated, in wider reaching but less 
readily interpreted fashion, by our cross-facility sensitivity com- 
parison (Fig. 5). 

In addition to carrying out the real-time multi-facility transient 
searches needed to realize these and other science gains, and en- 
abling follow-up EM observations by distributing transient alerts 
to interested observers, AMON will provide a framework for large 
observatory collaborations to work together on targeted archival 
analyses using AMON’s comprehensive events database and coinci- 
dence analysis toolkit. In these several ways, AMON will leverage 
and enhance the capabilities of existing and future high-energy 
and multimessenger observatories, powering a new and ambitious 
exploration of the transient universe using all four forces, and so 
helping to realize the immense promise of this dawning age of 
multimessenger astronomy. 
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Appendix A. Effective volume integrals 

In a number of examples in the text, a simple population model 
was applied that assumes a uniform spatial distribution of sources 
with a fiducial set of source parameters (encapsulated here by the 
parameter 2). In the case where z -c 1, we can define an effective 


volume V (eff) for a specific set of trigger conditions by the Euclidean 
volume integral, 


dV (eff) 

dQ 


1 f D 

-- — j 47zr 2 drP(true|/L r), 


(A.l) 


where P(truejr, A) is the conditional probability of a true positive 
detection, given the fiducial source parameters and distance to 
the source r. The maximum distance D may represent either a phys- 
ical limitation of the method or else a parameter that we use to reg- 
ulate the integral before taking D — > oo. The parameter(s) A may 
depend on pointing direction, and so we have written V <eff| as differ- 
ential with respect to solid angle Q. The expected rate of source 
detections is then given, after integrating over solid angle, by 

R = pV m , (A. 2) 


where p is the intrinsic rate per unit volume of transient sources. In 
the special case where the observatory is equally sensitive across its 
field of view, is independent of pointing and we can write 


R = pQ 


dV m 
da ' 


(A.3) 


In the case where the participating observatories are each mon- 
itoring for Poisson processes, P can be taken to be the cumulative 
Poisson probability distribution above some set of number thresh- 
olds n 0 . The vector notation includes the possibility of studying 
events from multiple observatories, with indices a = 1 -+ M, so 
that 


M f n 0a 

P(n > nofir) = n 1 - 

a = 1 \ k a = 0 


k n \ 


(A.4) 


where A a = i J dE^-A ( a elT> (E) is a measure of the number of parti- 
cles that can be detected by observatory a (where N a is the total 
number of particles produced and A[, eff| is the observatory’s effective 
area). It is normalized in such a way that 2 a /r 2 is the expected num- 
ber of detections for a source at a distance r. The trigger condition 
n n 0 is shorthand for simultaneous detections above the specified 
number threshold for each observatory. 

As such, we are interested in integrals of the form 



(A.5) 


where we have changed the integration variable to x 
fined A = K = £J a k 0 , and the term Qj(l) — 

By the multinomial theorem, we note that 


= A/r 2 and de- 

(q^)n 0 (#- 


V K (A)= ]T ( A - 6 ) 

,= K 

where we have used a single subscript on the left hand side (rather 
than a vector). As expected, this is the same result as Eq. A.5 applied 
to a single observatory. The interpretation is that Eq. A.6 can be ap- 
plied to triggers from multiple observatories that are interchange- 
able (e.g., when considering k = 2 neutrinos form IceCube or 
ANTARES, regardless of which of the observatories triggered). 

Importantly, v 0 and V\ diverge as D — > oo, but the linear combi- 
nation ^ - v 0 - V\ is finite. Thus, if we are looking for a single Pois- 
son trigger, then a physical limit D must be applied to achieve a 
finite result (e.g., by requiring that a single detected neutrino orig- 


68 


M.W.E. Smith et al. / Astroparticle Physics 45 (2013) 56-70 


inate from within the sensing region of a gravitational wave net- 
work). However, if we instead require that k > 2, then the result 
is finite as D — > oo. 

We can also marginalize over one observatory, defining 
k = (ki , k') and using the normalization of the Poisson series to find 

OO 

= (A.7) 

If, =o 


We can apply this result to two observatories. In the case where we 
search for at least m events from the first observatory and at least n 
events from the second, 


dV {e!f) 

Llv mn 

dQ 


D 3 

T 


ij^m,n 


q 3 m - 1 n-1 m-ln-1 

■y ~ E^oo - E y jfe) + EE y (jj) (^i.7 2 )- 

i=0 j=0 i = 0 j = 0 


(A.8) 


However, if we search for a total of n events, regardless of which 
observatory they come from, 


The photon spectrum is obtained by convolving a flat photon 
energy distribution from 7t° decay with the 2 1/3 0.97- ^ (where 
factor 2 accounts for the number of photons produced from each 
decaying pion and 1/3 is a fraction of neutral pions to the total 
number of pions) and the Hawking primary spectra. 

From decaying and n 1 , v ; , and v,, are produced, together 
with /( and fi + . These neutrinos have flat energy distribution, 
but there are also secondary neutrinos originating from the consec- 
utive muon decays: /z+H — > e +( -) + v e (v e ) + Vjj(Vp). Energy distri- 
butions of each of these neutrinos are convolved with 2/3-0.97- 
d ^ x and the quark/gluon Hawking spectra, and summed afterwards. 
Since production of v T from decaying hadrons is highly suppressed, 
most of this neutrino flavor originate from the direct Hawking radi- 
ation and from the decays of the directly emitted r leptons. There- 
fore, one would expect a flux of v x that is almost two orders of 
magnitude lower than the total flux of other neutrino flavors. We 
neglected these two v T contributions to the total neutrino flux in 
this work. 

B.2. SUSY decay chains 


dvf f) _ D 3 
dQ 


n-1 


Ew 

J=0 


(A.9) 


Both Eqs. A.8 and A.9 lead to finite results as D — > oo. 

We may wish to consider the case where one or both observa- 
tories have an effective area that changes significantly off-axis, 
meaning that depends on direction. In general, one must carry 
out the integration over solid angle numerically, although we have 
already noted the trivial analytic result in the case of Eq. A.3. In 
addition, there is an analytic solution in the case where only the 
first observatory has a dependence on zenith angle 0 Z , that depen- 
dence is proportional to cos 0 Z , and the number threshold of the 
first observatory is n, > 1. However, the result is sufficiently com- 
plicated that we choose to omit it here. 


Appendix B. Primordial black hole source model 

Here we describe the main decay chains leading to the final par- 
ticle spectra from the PBH explosion under the assumptions of the 
Standard Model (SM) and Supersymmetry (SUSY). 


If the Minimal Supersymmetric Standard Model (MSSM) de- 
scribes Nature at the high-energy scales associated with the explo- 
sions of PBHs, the number of dof available to be radiated increases 
by more than a factor of two: 244 compared to 118 dof in the SM, 
where we neglected graviton and gravitino dof and included five 
physical Higgs fields (see, for example, Ref. [113]). This leads to 
an increase of the factor a (M) that directly affects the rate of evap- 
oration by a factor of ~3.3 when the PBH temperature reaches the 
SUSY particle production energy scale. From that point, the time 
left until complete evaporation will be shorter compared to the 
time predicted under the SM assumption. For example, one second 
before complete evaporation, the temperature of the MSSM PBH 
would be T ~ 5.5 TeV with almost 50% more mass to be radiated 
than in the case of the SM PBH with T ~ 8 TeV at the same time left 
before its complete evaporation. 

Given a number of unknown supersymmetric parameters (more 
than 100), for simplicity, we assumed that gluinos (g) are heavier 
than squarks (q) (as in, for example, the SPSla benchmark scenario 
[114] and the mSUGRA B benchmark model [115]), thus they decay 
into antisquark/quark and squark/antiquark pairs [116]: 


BA. SM decay chains 

During the final stage of the PBH life, all of the 118 SM particle 
degrees of freedom (dof) are radiated away. Since there are more 
dof for quark and gluons (72) than for leptons and photons (26), 
the final spectra of both HE and UHE gamma rays and neutrinos 
originate mostly from the decaying hadrons [104]. Neutrons and 
antineutrons are produced from these decaying quark fragmenta- 
tion products, as well. 

We assumed that quark and gluon jets emitted during the final 
stage of PBH evaporation fragment into pions (7t°,7tr and 7t+) and 
baryons ( p.p , n, and n) with branching ratios of 0.97 and 0.03 [112], 
respectively. A fragmentation function d ^ x of quarks or gluons into 
a particle of type X has a form [112]: 

^ = H*-' 5 ( l-*) 2 , (A.10) 

where x = E X /E,E X is the energy of a particle of type X, and E is the 
quark/gluon jet energy. 

For baryon spectra, we assumed that equal numbers of p.p.n, 
and n are produced after hadronization. Thus, the final n+h spec- 
trum is obtained by multiplying distribution ^ by 0.015, convolv- 
ing it with the Hawking primary spectra for quarks and gluons and 
integrating over time. 


g^q + q.q + q. (A.ll) 

The right chiral states of squarks would decay then mostly into the 
lightest netralino (Z?): 

+ (A.12) 

whereas the left handed states would decay into charginos (%f) or 
heavier neutralinos (Z°) [117]: 

+ Z° + q- (A- 13) 

Further down these decay chains, we assumed that and jjjj pref- 
erably decay into leptons and charged leptons plus the lightest neu- 
tralino, respectively: 

zi - I' V; X °2 - l + rx°v (A.14) 

Without carrying a complete calculation, our rough estimate for the 
multimessenger particle fluxes under the assumption of the MSSM 
are: ~4 of the SM flux for neutrinos and ~3.5 of the SM flux for pho- 
tons and neutrons/antineutrons. 
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